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^ ■ Abstract. The self- localized quasi-particle excitation of the electron-positron field (EPF) is 

/""\ , found for the first time in the framework of a standard form of the quantum electrodynamics. 

This state is interpreted as the "physical" electron (positron) and it allows one to solve the 
t"*** \ following problems: i) to express the "primary" charge eo and the mass mo of the "bare" 

electron in terms of the observed values of e and m of the "physical" electron without 
, any infinite parameters and by essentially nonperturbative way; ii) to consider /x-mcson as 

Q-c another self-localized EPF state and to estimate the ratio m M /m; iii) to prove that the self- 

^ (— i localized state is Lorentz-invariant and its energy spectrum corresponds to the relativistic 

free particle with the observed mass m; iv) to show that the expansion in a power of the 
observed charge e <C 1 corresponds to the strong coupling expansion in a power of the 
"primary" charge e$ 1 ~ e when the interaction between the "physical" electron and the 
transverse electromagnetic field is considered by means of the perturbation theory and all 
t— I ■ terms of this series are free from the ultraviolet divergence. 

>' 
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!>■ ! 1 Introduction 

O ■ 

It is no doubt at present that the Standard Model is the fundamental basis for the theory of the 
electro- weak interaction pQ. It means that the quantum electrodynamics (QED) is actually the 
part of the general gauge theory. Nevertheless, QED considered by itself as the isolated system 
remains the most successful quantum field model that allows one to calculate the observed 
characteristics of the electromagnetic processes with a unique accuracy (for example, [3]). 
It is well known that these calculations are based on the series of rules connected with the 
perturbation theory in the observed charge e of the "physical" electron and the renormalization 
property of QED. The latter one means that the "primary" parameters of the theory (the charge 
eo and the mass too of the "bare" electron) that are defined by the divergent integrals, can be 
excluded from the observed values. However, even the creators of the present form of QED 
were not satisfied because "the calculation rules of QED are badly adjusted with the logical 
foundations of quantum mechanics and they cannot be considered as the satisfactory solution 
of the difficulties" [H § 81] and "it is simply a way to sweep the difficulties under the rug" [5]. 

There are a number motivations for calculation of the "bare" electron characteristics eo 
and toq in spite the fact that these values are unobserved. First of all it is the question whether 
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the system of Maxwell and Dirac equations as the mathematical model for the quantum field 
system is closed and self-consistent when considering the processes related to the interaction of 
electrons, positrons and photons? In that sense QED in the existing form of renormalization 
is an unclosed theory because it includes an additional, external dimensional parameter which 
enables the regularization of integrals. The contradiction between the small value of the coupling 
constant and its infinite calculated value shows a logical inconsistency of the perturbation theory 
in QED, so called "Landau pole" [6j § 128]. It is very important to consider the renormalization 
problem in the framework of the "logical principles of quantum mechanics" in order to develop 
some new nonperturbative approach in the field theory in application to the real physical system 
with completely defined Hamiltonian. At present nonperturbative methods are mainly studied 
for quite abstract quantum field models with a strong interaction (for example, [7]). Such 
methods may be especially interesting for the non-renormalized quantum field models. 

It is also very essential that the dynamical description of the internal structure of the "phys- 
ical" electron gives the fundamental possibility to consider /i-meson as an excited state of the 
electron-positron field as it has been shown by Dirac [8]. 

The relation between the "primary" coupling constant eo and the charge e is undetermined 
in the present form of QED. Therefore it is possible that the value eo is large in spite the ob- 
served renormalized charge being small eCl. Precisely this possibility (eo > 1, but e 1) is 
investigated in the present paper in order to find a spectrum of the quasi-particle excitations in 
QED without the perturbation theory. Our main goal is to find such a form of the renormaliza- 
tion that would be logically consistent but the calculation possibilities of QED for the observed 
values would be preserved. 

It is important to stress that the canonical QED model is considered with a nonzero mass of 
the "bare" electron too ^ and without the chiral symmetry. This approach is distinguished 
essentially from the nonperturbative analysis of the strong coupling QED model with a zero 
mass of the fermion field (for example, [9j [10] and citation therein). In the latter case the 
observed particle mass appears as the result of spontaneous breaking of the chiral symmetry, 
however, the cut-off momentum L is another undefined parameter in such theories. 

The article is structured in the following way. In Section [2] it is shown for the first time that 
the self-localized one-particle excitation can be found in the spectrum of QED Hamiltonian. 
This state cannot be calculated by means of the perturbation theory as a power series of the 
coupling constant eo- The stability of this state is conditioned by a localized charge distribution 
of the electron-positron field coupled with a scalar component of the electromagnetic field. 
The system of nonlinear equations for these spatially localized distributions is derived and its 
numerical solution is obtained. 

In Section [3] the self-localized excitation is interpreted as the "physical" electron with the 
observed values of the charge e and mass to. It allowed us to express characteristics eo and mo 
of the "bare" electron that are actually the parameters of the initial Hamiltonian in terms of e 
and m. It is shown that the relation between these values includes the singularity in the limit 
of e — > and cannot be calculated by means of the perturbation theory. This result cannot be 
also obtained in the framework of the "quenched QED" model based on the Schwinger-Dyson 
equation [11] or on the variational approach [12] because the charge renormalization did not 
take into account in this model. 

The considered physical interpretation of the self-localized state leads also to an important 
consequence. It is shown in Section [3] that there is another one-particle excitation of the electron- 
positron field with the same charge and spin as for electron but with the larger mass. Following 
Dirac [5] this excitation was considered as the "physical" //-meson and the ratio m/m^ is cal- 
culated. The calculated mass of /i-meson proved to be very close to its experimental value. It 
is essential that unlike [8] the //-meson mass is calculated without any additional parameters of 
the model with the exception of eo and toq. 
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The localized charge distribution in the "physical" electron corresponds to the spontaneous 
Lorentz symmetry breaking for the considered system. This phenomenon is typical for the 
quantum field theories with the particle-field strong coupling (for example, the "polaron" prob- 
lem \13\ I14j). It is shown in Section [5] that reconstruction of this symmetry leads to the depen- 
dence of energy of the excitation on its total momentum. It is essential that this dependence 
corresponds exactly the relativistic kinematics of a free particle with the observed mass m. 

It is shown in Section [5] that the strong coupling series in the initial QED Hamiltonian cor- 
responds to the perturbation theory in terms of observed charge e ~ eg < 1, However, all 
high-order corrections of this perturbation theory is defined by the convergent integrals with- 
out any additional cut-off parameters. The correspondence of these results with the standard 
renormalization procedure when calculating the observed characteristics of the electromagnetic 
processes is also discussed. 

2 Self-localized state with zero momentum 

It is well known that the spatially localized states are very important for a lot of quantum field 
models. Let us remind the one-dimensional model for the scalar field with the Hamiltonian (for 
example, [15]): 

*-H{HT)'M«r-?)}- 

with the field operators <~p{x) and the corresponding momentum operators tt. 

This Hamiltonian has the eigenvector \tp) with the energy density localized in the vicinity 
of an arbitrary point xq. With the zero total momentum of the system the energy density is 
distributed as follows: 

(1) 

and corresponds to the finite total energy E = (2>/2/3) (m 3 /A). Lorentz -invariance of the 
system leads to the relativistic dependence of the energy on the momentum of the localized 
excitation [15j . 

It is important for the further discussion that the state (JTJ) cannot be derived by means of the 
perturbation theory based on the coupling constant A connected with the nonlinear interaction. 
One should separate the non zero classical component from the field operators y>(x) = {(p\ip\(p) 
averaged over the considered state even in the zeroth-order approximation. And the classical 
function ip(x) is not reduced to the plane wave as in the zeroth-order perturbation theory. It is 
satisfied to the nonlinear differential equation that is defined by the variation of the classical 
functional. 

The Frdhlich model for "polaron" problem gives another example |16j . This model corre- 
sponds to the electron-phonon interaction in the ionic crystals and is described by the following 
Hamiltonian: 

» = ^+E 4«j + 23/4 \/f £ fa + 

k k 

In this case the spatially localized electron state ("polaron") cannot be found also by means 
of the perturbation theory in terms of the electron-phonon coupling constant a. The variational 
wave function of the electron and the classical part of the phonon field = (a?) should be 
calculated even in the zeroth-order approximation. The translational symmetry of the initial 
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Hamiltonian leads to the dependence of the energy of this one-particle excitation on its total 
momentum |14j . 

Let us now consider the nonperturbative analysis of the spectrum of the one-particle excita- 
tions of the QED Hamiltonian that is defined by the following form (for example, [17]): 

H = j dr:{ip*(r)[a(p + e A(r)) +/3m ]^(r) + e <p(f)p(r) - ^(V(^(r)) 2 }: +^uj(k)n^ x , 

kX 

m = \ w (f)^(f) - (oi ■ ( 2 ) 

We suppose here that the field operators are given in the Schrodinger representation, the 
spinor components of the electron-positron operators being defined in the standard way [T7] 

In these formulas h = c = 1; the primary charge (— eo), &o > and mo are considered as the 
parameters of the model; the symbol : H : means the normal ordering of the operators excluding 
the vacuum energy [3]; a, (5 are Dirac matrices; u^ sv and v^ su are the components of the bispinors 
corresponding to the solutions of Dirac equation for the free "bare" electron and positron with 
the momentum p and spin s; ap s (a^ s ) and £y s (6^t ) are the annihilation (creation) operators for 

the "bare" electrons and positrons in the corresponding states. The field operator A(r) and the 
operator of the photon number n^ A are related to the transversal electromagnetic field and their 
explicit form will be written below. 

This Hamiltonian corresponds to the Coulomb gauge [T7], when the electron-positron field 
interacts with the scalar field and with the transversal photons of the electromagnetic field. It 
will be shown below that the reconstruction of the QED gauge symmetry connected with the 
longitudinal field does not change the form of the considered one-particle excitation. In the 
Coulomb gauge the operators of the scalar field 

0(f) = V4n J dk0p^ F 

can be excluded from the Hamiltonian [T7] . For that purpose one should use the solution of the 
operator equations of motion for ipr assuming that the "bare" electrons are point-like particles 
and "self-action" is equivalent to the substitution of the initial mass for the renormalized one. As 
a result the terms with scalar fields in the Hamiltonian are reduced to the Coulomb interaction 
between the charged particles. However, this transformation of the Hamiltonian ([2]) can not be 
used in this paper because only the dynamics of the mass renormalization is the subject under 
investigation. 

There is another problem connected with a negative sign of the term corresponding to the 
self-energy of the scalar field. If the non-relativistic problems were considered then the operator 
of the particle kinetic energy would be positively defined and the negative operator with the 
square-law dependence on <p(r) would lead to the "fall on the center" [18] as the energy minimum 
would be reached at an infinitely large field amplitude. However, if the relativistic fermion field 
is considered then the operator of the free particle energy (the first term in formula ([2])) is 
not positively defined. Besides, the states of the system with the negative energy are filled. 
Therefore, the stable state of the system corresponds to the energy extremum(!) (the minimum 
one for electron and the maximum one for positron excited states). It can be reached at the 
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finite value of the field amplitude (see below). The same reasons enable one to successfully 
use the states with indefinite metric [2] in QED although it leads to some difficulties in the 
non-relativistic quantum mechanics. 

According to our main assumption about the large value of the initial coupling constant eo 
we are to realize the nonperturbative description of the excited state which is the nearest to the 
vacuum state of the system. The basic method for the nonperturbative estimation of the energy 
is the variational approach with some trial state vector for the approximate description 
of the one-particle excitation. The qualitative properties of the self-consistent excitation in the 
strong coupling limit [13] show that such trial vector should correspond to the general form 
of the wave packet formed by the one-particle excitations of the "bare" electron-positron field. 
Besides, the effect of polarization and the appearance of the electrostatic field <p(f) should be 
taken into account, so we consider |<l>i) to be the eigenvector for the operator of the scalar field. 
Now, let us introduce the following trial vector depending on the set of variational classical 
functions U^ s , V^ s , </?(r) for the approximate description of the quasi-particle excited state of the 
system: 

|*i> ~ |$f (U&VfrM?))) = j dq{U<tsa± +^6+}|0;0;^(r)), 

£(f)|0; 0; <p(r)) = <p(?)\0; 0; ^(f)>, o r ,|0; 0; <p(r)) = 6 r ,|0; 0; <p{?)) = 0. (3) 

The ground state of the system is |<J>o) — 1 0; 0; 0) , if we use the same notation. It corresponds 
to the vacuum of both interacting fields. 

Firstly, let us consider the excitation with the zero total momentum. Then the constructed 
trial vector should satisfy the normalized conditions resulting from the definition of the total 
momentum P and the observed charge e of the "physical" particle: 

J2dQ[\U qs \ 2 + \V qs \ 2 } = 1, (4) 

s 

(5) 

The condition requires that the functions U qs and V qs should depend on the modulus of the 
vector q only. Besides, one should take into account that the trial vector } is not the accurate 

eigenvector of the exact integrals of motion Q and P as it represents the accurate eigenvector of 
the Hamiltonian |<3?i) only approximately. Therefore, in the considered zero approximation the 
conservation laws for momentum and charge can be satisfied only on average, and this leads to 
the above written normalized conditions. Generally, equation ([5]) should not be considered as the 
additional condition for the variational parameters but as the definition of the observed charge 
of the "physical" particle at the given value of the initial charge of the "bare" particle. Therefore 
the sign of the observed charge is not fixed a priori. Calculating the sequential approximations 
to the exact state vector (see Section H|) should restore the accurate integral of motion as 
well. An analogous problem appears in the "polaron" theory when the momentum conservation 
law was taken into account for the case of the strong coupling (for example, [14^ I19j). 

The trial vector is actually the collective excitation of the system and in this respect the 
variational approach differs greatly from the perturbation theory where the zero approximation 
for a one-particle state is described by one of the following state vectors: 

lM PT) e) = 4J0;0;0), ~ \<S>[ PT) p) = &jj0;0;0). (6) 

These vectors do not depend on any parameters and are eigenvectors of the momentum and 
charge operators. But they correspond to one-particle excitations determined by the charge eo of 



(d>( o )|P|d>f)) = j2<mw!s\ 2 + i^i 2 ] = P = o, 

s 

(<I>S 0) |Q|$S 0) ) = e "£dq[\V qs \ 2 - \U qs \ 2 ] = e. 

s 
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the "bare" electron and the field f(r) = 0. We suppose that the introduction of the variational 
parameters into the wave function of the zero approximation will enable us to take into account 
the vacuum polarization. 

It should be noticed that the another reason of inconsistency of the states © as the "physical" 
electron states because of the low frequency photon field was considered recently in [20] . 

So, the following variational estimation for the energy -Ei(O) = E\ (P = 0) of the state 
corresponding to the "physical" quasi-particle excitation of the whole system is considered in 
the strong coupling zero approximation: 

£x(0) * E[ 0) [U qs ;V qs ;cp(r)} = (^\H\^\ (7) 

where the average is calculated with the full Hamiltonian ([2]) and the functions U qs , V qs , (p(r) 
are to be found as the solutions of variational equations 

dEf\u qs -V qs -y(r)) dE{ 0) dE[ 0) 

dU qs dV qs dy{r) 1 ] 

with the additional conditions (|3J), (J3J). 

It is quite natural, that the ground state energy is calculated in the framework of the consid- 
ered approximation as follows: 

E ~ E^ = ($ \H\<1> ) = 0. 

Further discussion is needed in connection with the application of the variational princip- 
le (0, ([8]) for estimating the energy of the excited state, because usually the variational principle 
is used for estimating the ground state energy only. As far as we know it was first applied 
in [21] for nonperturbative calculation of the excited states of the anharmonic oscillator with 
an arbitrary value of anharmonicity. This approach was called the "principle of the minimal 
sensitivity" . It was shown in our paper [22] , that the application of the variational principle to 
the excited states is actually the consequence of the fact that the set of eigenvalues for the full 
Hamiltonian does not depend on the choice of the representation for eigenfunctions. As a result, 
the operator method for solving Schrodinger equation was developed as the regular procedure for 
calculating further corrections to the zero-order approximation. Later this method was applied 
to a number of real physical systems and proved to be very effective when calculating the energy 
spectrum in the wide range of the Hamiltonian parameters and quantum numbers ([23 1 1241 125j 
and the cited references). 

The average value in equation ([7|) is calculated neglecting the classical components of the 
vector field. They appear in the high-order corrections that are defined by the renormalized 
charge e<l and can be considered by means of the canonical perturbation theory (Section [5]). 
It means that 

($(°)|^*(r)[ai'(r)]V»(r)|^S 0) ) = 0. 

It should be noted that the possibility of constructing self-consistently the renormalized 
QED at the non-zero vacuum value of the scalar field operator was considered before [26] but 
the solution of the corresponding equations was not discussed. 

Then the functional for determining the zero approximation for the energy of the one-particle 
excitation is defined as follows: 



^ ' ^ ' s,s' H,v 

- tW«fv M [(ag + 0mo)„u + e^S^V^v^e^-^ - \ f df[V^(r)} 2 . (9) 
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In order to vary the introduced functional let us define the spinor wave functions (not opera- 
tors) which describe the coordinate representation for the electron and positron wave packets in 
the state vector |$^° ): 

*" ^ = / j&fa ? u ^ e - w = I wf E y>^. (io) 

In particular, if the trial state vector is chosen in one of the forms ([6]) of the standard 
perturbation theory, the wave functions (|10p coincide with the plane wave solutions of the free 
Dirac equation. For a general case the variation of the functional ([9]) by the scalar field leads to 



E(0) = I dr{^ + (r) 
^ +c (r) 



-iaV + /3m ) + -e <p(r) 



I 



(-iaV + 0m o ) + ^e (p(f) 



*(r) 
,C (r), 



dr[y + (r)^(r) + ^ +c (r')^ c (f ')] = 1, (11) 



ip{f) = £° / [$ + (?')y(f') - * +c (r')^ c (f')]. (12) 



4tt J \r — f' : 

The main condition for the existence of the considered nonperturbative excitation in QED 
is defined by the extremum of the functional (jlip corresponding to a non-zero classical field. 
The structure of this functional shows that such solutions of the variational equations could 
appear only if the trial state vector simultaneously included the superposition of the electron 
and positron wave packets. So, such solutions cannot be obtained by means of the perturbation 
theory with the state vectors j6|). 

Equation (|llj) and the Fourier representation ([3]) for the trial vector clearly indicate that the 
assumption concerning the localization of the functions \&(r) near some point does not contradict 
to the translational symmetry of the system because this point by itself can be situated at any 
point of the full space with equal probability. The general analysis of the correlation between 
the local breaking of the symmetry and the conservation of accurate integral of motion for the 
arbitrary quantum system was considered in detail by Bogoluibov in his widely known paper 
"On quasi-averages" [27]. A similar analysis of the problem in question will be given in Section [H 

Varying the functional (llip by the wave functions ^(r) and x I /C (r) taking into account 
their normalization conditions leads to the equivalent Dirac equations describing the electron 
(positron) motion in the field of potential ip(r): 

{(-iaV + [3m ) + e ip(r)}y(r) = 0, {(-iaV + (3m ) + e ^(f)}* c (r) = 0. (13) 

But it is important that in spite of the normalization condition (llip for the total state 
vector (jlOp each of its components could be normalized differently 



/ 



dr^(r)^(r) = — J dr^ c (r')^ c (r') = (14) 



The constant C is an arbitrary value up to now. It defines the ratio of two charge states 
in the considered wave packet. As a result the self-consistent potential (f(r) of the scalar field 
depends on C because of the equation (fT2l . 

We should discuss the procedure of separating variables in more detail, because of the non- 
linearity of the obtained system of equations for the wave functions and the self-consistent 
potential. Since the considered physical system has no preferred vectors if P = 0, it is natural 
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to regard the self-consistent potential as spherically symmetrical. Then the variable separation 
for the Dirac equation is realized on the basis of the well known spherical bispinors [2] : 



jlM 



g{r)SljiM 
if(r)£lji> M 

Here Q,um ar e the spherical spinors [2] describing the spin and angular dependence of the 
one-particle excitation wave functions; j, M are the total excitation momentum and its projec- 
tion respectively, the orbital momentum eigenvalues are connected by the correlation l + l' = 2j. 
It is natural to consider the state with the minimal energy as the most symmetrical one, cor- 
responding to the values j = 1/2, M = ±1/2, I = 0; I' = 1. This choice corresponds to the 
condition according to which in the non-relativistic limit the "large" component of the bispinor 
*I> ~ g corresponds to the electronic zone of the electron-positron field. Then the unknown 
radial functions /, g satisfy the following system of the equations: 

^j^- ~ -( r 9) ~ ( m o ~ e <p{r))(rf) = 0, 
ar r 

&1 + I (r/ ) _ ( mo + eo<p(r))(rg) = 0. (15) 
ar r 

The states with various projections of the total momentum should be equally populated 
in order to be consistent with the assumption of the potential spherical symmetry with the 
equation (112p . So, the total wave function of the "electronic" component of the quasi-particlc 
excitation of the electron positron field is chosen in the following form: 



* = 7f[*l/W/2 + *l/2A-l/ 2 ] = ( ) 



X, + = ^[n l /2 |Jl l/2+ n i/2.i ) -l/2], 1 = 0,1- (16) 

In its turn, the wave function \& c is defined on the basis of the following bispinor: 

•*» - ( 7;<Sr ) • « i7 » 

The radial wave functions fi, g\ in this case satisfy the following system of equations 

^EA + \ rgi) _ ( mo + eoV (r)){rh) = 0, 
ar r 

<MA _ i( r/l ) _ (mo _ eo <p(r))(r gi ) = 0. (18) 
ar r 

These equations correspond to the positronic branch of the electron-positron field with the 
"large" component ~ g\ in the non-relativistic limit. 

It is important to note that the functions and ^ c satisfy the equations (|13|) that have the 
same form. It imposes an additional condition of the orthogonality on them: 

(^ c |$)=0. (19) 

Taking into account this condition and also the requirement that the states with different 
values of M should be equally populated one finds the "positronic" wave function 



71^/2,0,1/2- *l/2 ,0,-l/2\-[ gi{r)x - 



Xi =^[^1/2,2,1/2-^1/2^-1/2]) 1 = 0,1. (20) 
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The equation for the self-consistent potential follows from the definition of <p(r) in for- 
mula (|12|) taking into account the normalization of the spherical spinors [2J: 

(ftp 2 dip e r 2 2 2 2i 

The boundary condition for the potential is equivalent to the normalization condition (J5J) and 
defines the charge e of the "physical" electron (positron) 

^(r)|™ = J- = £L I™ rldnifir^+g^n) - /f (n) - gUn)]. (22) 

It is important to stress that the form of the functions given above is defined practically 
uniquely by the imposed conditions. At the same time the obtained equations are consistent 
with the symmetries defined by the physical properties of the system. The first symmetry is 
quite evident and relates to the fact that the excitation energy does not depend on the choice 
of the quantization axis of the total angular momentum. 

Moreover, these equations satisfy the condition of the charge symmetry [2j. Indeed, by direct 
substitution, one can check that one more pair of bispinors leads to the equations completely 
coinciding with (fl5j) . (fP8l) 



*jlM=\ ), (23) 



It means that these bispinors allow one to find another pair of the wave functions which are 
orthogonal to each other and to the functions (|16p . (|20p but include the same set of the radial 
functions 

a / Wi(r)Xi \ ^c_( ~f(r)xt \ (ok) 
^ " V -AWxo ) ' V -W(r)xt ) ■ (25) 



These functions differ from the set (|16|) . (j20|) because they lead to a different sign of the 
observed charge of the quasi-particle due to the boundary condition fl22f) and describe the "physi- 
cal" positron. 

The structure of the equation f)13[> shows that the considered variational method is consistent 
with the gauge symmetry of the initial Hamiltonian. One can see that these equations are 
invariant with respect to the following transformations: 

V =^ V + ie Ai(r), 

e _i/J W*(f), V/3(f) = eoAiir), 

with an arbitrary longitudinal potential Ai(r). 

It means that the Hamiltonian ([2]) could be chosen in an arbitrary Lorentz gauge with the 
classical components both for the scalar field <p(f) and for the longitudinal field Ai(f) if the 
following condition was fulfilled: 

A/3(f) = pi(r) = 0, pi(r) = ^ + (f)(f - V)*(r). 

One can easily check that the condition pi (r ) = is fulfilled identically for the functions that 
describes the quasi-particle above. 

Let us now proceed to the solution of the variational equations. It follows from the qualitative 
analysis that the important property of the trial state vector is the possibility to vary the relative 
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contribution of the electronic and positronic components of the wave function. Therefore let us 
introduce the variational parameter C in the following way: 



r 2 dr[f 2 (r)+g 2 (r)} 



1 



o 



l + C 



r 2 dr[f 2 (r)+g 2 (r)] 



C 



o 



l + C 



The dimensionless variables and new functions can be introduced 



x = rrriQ, E = eniQ, eo^(r) = mo4>(x), 
v{x)y/m^ = rf(r), u 1 (x) y /mo' = rg x (r), 



a , 



4tt 

v 1 (x)^/m^ = r/i(r) 



u(x)y/mo = rg(r), 



(26) 



As a result the system of equations for describing the radial wave functions of the one- 
particle excitation of the electron-positron field and the self-consistent potential of the vacuum 
polarization can be obtained: 



du 1 

u — (1 — <p{x))v = 0, 

ax x 

du\ 1 , , . . . 

-^ + -u 1 -(l + tf>(x))v 1 = 0, 



dv 1 . , , 
— + —v - (1 + 4>{x))u 
dx x 



0, 



-j— ~ -vi - ( 1 - 4>(x))ui 
dx x 



<f>(x) = a 



y 



X 



dyp{y) 



0, (27) 
p(x) = [u 2 (x) + v 2 (x) — u\(x) — v\(x)\ . 



The mathematical structure of equations (|27p is analogous to that of the self-consistent 
equations for localized state of "polaron" in the strong coupling limit [13] • Therefore the same 
approach can be used for the numerical solution of these nonlinear equations. It has been 
developed and applied [28] for the "polaron" problem on the basis of the continuous analog of 
Newton's method. 

Let us take into account that the system of equations (127j) can be simplified because the 
pairs of the functions u, v and u\, V\ are satisfied by the same equations and differ only by the 
normalized condition. Therefore they can be represented by means of one pair of functions if 
special notations are used: 



u 



1 



;u , 



l + C 

dx[ul(x) + Vq(x)] 
duo 1 



1 



-,v , 



C 



i + c u ' V l + c 

1, p (x) = Uq{x) +Vq(x), 



vo, 



Vl 



C 



l + C 



dx 
</)(x) 



-u - (1 - <f>(x))v = 0, 



dv , 1 

dx x 



oco 



l-C 
l + C 



4>o(x), 4>o(x) 



dy 



(1 + 0(x))«o = 0, 

Po(y) , l f x , f , 

h - / dyp {y). 

y x J 



(28) 



The energy of the system (jlip can also be calculated by these functions: 



£i(0) = E(0) = m 



1 - C 
l + C 



m 1 1 - c 
T H — an — 

2 l + C 



n 



T 



dx[(u v - v' u ) 



t U Vo 



+ (v 2 o-v 2 )}, 



n 



dx(f> Q (uQ + Vq) 



(29) 



and equation (|28p can be obtained when varying of the functional (|29p . 

The required solutions are to be normalized and this condition defines the asymptotic be- 
havior of the functions near the integration interval boundaries: 



u 



Ax 



1 + 



1 



! (0) 



.4 



1 - 0(0) 



0. 



no ~ A\ exp(— x), 



vo 



-A\ exp(— x), 



oo. 
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Figure 1. Localized wave functions of the quasi- 
particle excitation; x = r/ro = 2 ^°^ r , u = 
&g(xr ), v = u = £xf(xro), £ = i(j^b) 3/2 - 



Figure 2. Self-consistent potential of the exci- 
tation; x = r/r = ^t^ 1 , <j>o(x) = ira^O^o)- 



The value 

1 -C 



is the free parameter of the equations (|28p and it plays a role of the eigenvalue when the nontrivial 
normalized solution exists. 

The method for the numerical solution of the nonlinear self-consistent system of the equa- 
tions (|28h was described in detail in the paper [29]. Only the numerical results for the localized 
wave functions and for the scalar potential are described in the present work. Certainly, the nu- 
merical value for the parameter a depends on the accuracy of the finite-difference approximation 
for the differential operators in the whole interval of integration. The value a was calculated 
more accurately in comparison with [29j : 



a = ao rs —3.531. (31) 

Fig. Q] shows the solutions uq, vq for the electron and positron components of the excitation 
that are localized in the domain with the linear size of ~ TOq . Fig. [2] represents the self- 
consistent potential (fro that provides stability of the system in this domain and corresponds to 
the value do- It gets over the Coulomb potential of the "physical" charge e for r > ro = rriQ 1 . 
It is important that the characteristic size of this excitation ro is the same order as the classical 
radius of the electron r e = a/m, namely ro = gfen ~ 0.15r e (see below equation (|38p). 



3 Physical interpretation of the quasi-particle excitation 
and estimation of the /x-meson mass 

The stationary localized collective excitation of the electron-positron field described above is 
of great interest by itself as the eigenvector of the well known QED Hamiltonian that cannot 
be calculated by means of the perturbation theory and has not be considered before. But it 
is also essential to find its physical interpretation because the only stable objects observed in 
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the electrodynamic processes are electrons (positrons) and photons. Therefore it is natural to 
suppose that this localized state describes the "physical" electron (positron) with the observed 
charge e. The integral charge of the considered one-particle excitation is defined by the boundary 
condition (|22p and this supposition leads to: 



e »irr^H (32) 

Taking into account the definition (|3U|) one can found the following relation between the 
"primary" coupling constant ao = e o/4vr and the observed value of the fine structure constant 
a = e 2 /47r 

a 2 

a = ss 1708.1. (33) 
a 

This formula defines the renormalization of the charge in the considered approximation 
and shows self-consistency of the initial supposition that the interaction between the "pri- 
mary" electron-positron and scalar fields is strong. Then the renormalization constant [2] 
(a = Z(a)oto) is: 

= £ (34) 

It should be stressed once more that the large value of the "primary" coupling constant 
ao does not mean at all that the perturbation theory cannot be applied for the calculation of 
the observed physical values. In connection with it let us remind that the "primary" coupling 
constant in the existed form of QED tends to infinity (ao — ► oo) because of "Landau pole" [BJ 
§ 128]. Nevertheless, it can be excluded from the observed values by means of the renormalization 
procedure. In our representation the large but finite value of ao is important only for the 
formation of the initial basis of the self-localized states but it is also excluded when calculating 
the observed physical values in a power series of the "physical" coupling constant a ~ 1 <C 1. 
However, in this case one can avoid the divergent integrals when performing the renormalization 
procedure (see below Section EJ). 

The QED Hamiltonian (J2J) is defined with eo > and the integral charge of the quasi- 
particle e < because of the conditions f|30j) and (|3ip . It means that the considered excitation 
corresponds to the "physical" electron. The excitation corresponding to the "physical" positron 
is defined by the bispinors (f23|) . 

The change of the integral charge of the excitation is explained by the local intersection 
of the electron and positron energy zones in the strong self-consistent scalar electromagnetic 
field. The analogous states appear also in the strong Coulomb field of the nucleus with the 
large charge Z > 137 when the gap between Dirac zones tends to zero and stabilization of the 
electron states is achieved due to creation of the additional positrons (see, for example, [30| 
and references therein). In our case the origin of the charge density of one sign is the source of 
the strong scalar electromagnetic field that leads to coming together the electron and positron 
levels. This field can be compensated by the extra value of the charge density of opposite sign. 
More detailed qualitative analysis of the structure of the quasi-particle excitation was discussed 
in the paper [31]. It is also interesting to note that the relation between the charges of the 
"bare" and "physical" electrons has the same form as the equation for the electric and magnetic 
charges in the theory of the Dirac monopole [32]. 

The eigenvalue ao corresponds to minimum of the functional (|29p with fixed functions uq, 
vq, (fo with respect to the parameter £ = (1 — C)/(l + C) and can be expressed in terms of the 
integrals T and II that define kinetic and potential contributions in the total excitation energy 

T 

ao = -go (35) 
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and the numerical value of the integral T can be calculated with the considered accuracy as 

T « 0.749. (36) 

Then the total energy of the excitation with zero momentum is: 

. mn Tan T 
E(0) = °- = -m a— > 0. (37) 

OiQ Z ZCLq 

This value defines the minimal energy of the one-particle excitation of the electron-positron 
field and its positive sign corresponds to the "bottom" of the "physical" electron zone in the 
renormalized QED. It is also consistent with the negative charge of this excitation |31| . It is 
natural to consider this value as the rest-mass m of the "physical" electron and to find the mass 
renormalization in QED: 

T 2|oo| 

£■(0) = m = m e = —moa- — , tuq = m e ps 1291. 7m e . (38) 

This relation (E(0) <C mo ) shows that the primary mass of the "bare" electrons and 
positron ~ mj is compensated by the binding energy of their charge distributions almost 
completely. It is interesting that the characteristic size of the domain, where the considered 
one-particle excitation is localized (Ar ~ ), is the same order as the value of the electron 
"radius" (r e = a/m) in the classical model of Abraham-Lorentz [8]. 

In accordance with the formulas (|33j) , (|38p characteristics of the charge spatial distribution in 
the "physical" electron depend on the observable QED parameters a and m only. Such internal 
structure of the electron does not contradict to the well known results (for example, [33]) that 
maintain that the observed characteristics of the electron cannot depend on any dimensional 
parameter with the exception of m (see below the discussion in Section [5]). Actually the local 
charge distribution in QED arises as well in the framework of the perturbation theory when 
considering the physical interpretation of the renormalization of charge (for example, [2]). 

As it was shown by Dirac [8], investigation of the "physical" electron with the distributed 
charge is of great interest because it gives the possibility to interpret the "physical" /u-meson as 
the excited state of such system. In order to describe the dynamics of such excitation Dirac in- 
troduced the hypothetical elastic parameter that is absent in QED. However, the variational ap- 
proach considered in the present paper allows one to analyze the one-particle excitation differed 
from the "physical" electron without inclusion of any additional parameters. 

Let us choose the trial state vector in the same form as it is given by formulas 0, 0, 
but with the wave functions that are orthogonal to the "physical" electron state vector. These 
functions are satisfied to the system of the self-consistent equations in the form fll 1H — (j 14|) . but 
with the nonzero eigenvalues E\ ^ E<i for both components of the "primary" electron-positron 
field: 

{(-iaV + Pmo) + e ^(r) - E^^f) = 0, 
{{-iaV + pm ) + e ^(r) - E 2 }%{r) = 0, 

J dr^(r)^,(r) = — L^, J dr^ c (r')%(r') = (39) 
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In this case the energy that defines the observed mass of the "physical" /i-meson can be 
calculated as 

E^P = 0)=m„ = Ex- E 2 . 

The parameters C, eo, mo in equations (f38|) should be the same as was defined by equa- 
tions (f32j) - ([38|) because the observed charges of the "physical" electron and /n-meson coincide. 
Besides, the same bispinors (|16p -(|24 p as for electron should be used when separating variables 
in equation (|39p . However, the radial functions for both components of the "primary" electron- 
positron field that form the "physical" /i-meson should correspond to the different eigenvalues. 

If the dimensionless variables (|26|) have been used the following system of equations is ob- 
tained: 

^ - - (1 + ei - <^(x)K = 0, ^ + ^ - (1 - ei + M x )) u n = 0, 

^ + \u XlL - (1 - e 2 + ^(z))uip = 0, ^ - ^ - (1 + e 2 - <^(z))tti M = 0, (40) 



y x ' 



P(x) = [ulix) + uj(z) - u? M (s) - . 



If the solution of the nonlinear problem (|40p with the eigenvalues e\p and the orthogonality 
and normalization conditions is found, the observed mass of //-meson is defined by the following 
formula: 



m (Ip - I 



UnV,, , o O \ 1 



^ = y dx WpVp - u'^Vf, - 2-^ + («J - «J) + -oo0 M («J + «J) j , 

where integral Jia is expressed by the functions ttin, v\n with the same formula as the integral In 
by the functions u^, v^; the parameter ao is defined by equation (J35l) . 

Equation (|37p for the observed mass of the electron can be written in the same form 

1-C 

where integral Iq is expressed by the functions uq, vq from equations ([27]) . 

Calculation of the excited states for the nonlinear equations (|40p with the self-consistent 
potential proved to be a quite difficult numerical problem. In the present paper this solution 
was calculated on the basis of the "adiabatic" approximation for the potential. It means that 
the potential was fixed in the same form as in equations (|27p : (f>a ~ 4>q. In this case the functions 
Uu, Vn coincide with uq, vo but the functions uin, ui« correspond to the first excited state in 
the Dirac equations (|27p . Then the relation (|32p and the definitions of the integrals I lead to 
the following formulas: 

1 ~ C 6 T T 

1 + 6 e 

and the observed electron-//-meson mass ratio can be calculated as: 

m e ~ 1 + C 21 p 2a ~ ' / M 

The calculated values for the integrals are the following I M « 0.95, /i^ « 0.25 and the mass 
ratio is estimated as 

^ « 194. 
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This number can be compared with the experimental value (m^/m e )E ~ 206 and the result 
obtained by Dirac [B]: {m^/m^D ~ 54. Thus, the physical interpretation of the QED one- 
particle excitations leads to a quite good estimation for the observed mass of /i-meson. It can 
be improved if the completely self-consistent solution is found for equations (|40p . 

It may seem that interpretation of the "physical" /U-meson as the excited state of the "bare" 
electron-positron and electromagnetic fields contradicts to the experimentally observed conser- 
vation of the muon (not electrical) charge in the electromagnetic processes. However, let us 
remind that this state is the collective(!) excitation of the whole system with the scalar elec- 
tromagnetic field (p^ not equal to the field cp for the "physical" electron self-localized state. 
In accordance with the formula ([3]) it means that the both states are defined by the coher- 
ent states corresponding to different field amplitudes. Therefore the transition between these 
states is defined by the overlapping integral between the corresponding coherent states of the 
electromagnetic field. If one uses the standard definition of such states [2] this integral can be 
represented in the following form 



exp 



-4tt / —\E(k)-E^k)\ 
J u k 



(41) 



Here E(k), E^(k) are the Fourier components of the corresponding electric field strengths, 
for example: 

= - 7^)3 / drfiptfe** 
With utilisation of the dimensionless variables (|26p the integral (|4ip can be written as 

roc 

I = e- aoA , A= / tdt[$(t) - %(t)} 2 , 
J o 

roc f'OO 

= / xdx(j)'(x)sm(tx), = / xdxcf/(x)sm(tx). (42) 

Jo Jo 

Parameter A ~ 1 is defined by the converged integrals because of the asymptotic behavior of 
both functions 4>(x), (f>u(x) when x — > oo and x — > 0. The overlapping integral (|42|) permits one 
to estimate the electromagnetic lifetime of //-meson as 

- mo 1 /" 1 ~ 10 300 [sec]. 

This value ensures the conservation of the muon charge in the electromagnetic processes be- 
cause it is essentially bigger than the total /x-meson lifetime conditioned by the weak interaction. 



4 Lorentz invariance of the self-localized state 
with nonzero momentum 

In the previous sections the resting quasi-particle with a non-trivial self-consistent charge dis- 
tribution, the finite energy E(0) and a zero total momentum P = was considered in the 
framework of a nonperturbative QED. 

The obtained solution allows one to imagine the internal structure of the resting "physical" 
electron (positron) as a strongly coupled state of charge distributions of the opposite sign. The 
large values of integral charges of these distributions compensate each other almost completely 
and their heavy masses are "absorbed" by the binding energy. Actually the energy ±E(0) 
defines the boundaries of the renormalized electron and positron zones resulting from the strong 
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polarization of the electron-positron field when the excitation appears. But this excitation could 
be interpreted as the "physical" electron (positron) if the sequence of the levels in every zone 
determined by the vector P / were described by the Lorentz invariant relativistic energy 
spectrum ±E(P) of real particles, that is 



E{P) = y/P 2 + E 2 (0) = y/P2+m?. (43) 

It is worth saying, that the problem of studying of the dynamics of the self-localized excitation 
should be solved for any system with a strong interaction between quantum fields in order to 
calculate its effective mass. For example, a similar problem for Pekar "polaron" [13] in the ionic 
crystal was considered in \14:\ \19\ [34"1 [35] and in a lot of more recent works. It is essential that 
because of the non-linear coupling between the particle and a self-consistent field the energy 
dispersion E(P) for the quasi-particle proves to be very complicated. As the result, its dynamics 
in the crystal is similar to the motion of the point "physical" particle only at a small enough 
total momentum. 

However, in the case of QED the problem is formulated in a fundamentally different way. 
There is currently no doubt that the dynamics of the "physical" excitation should be described 
by the formula (|43p for any(!) values of the momentum P because of the Lorentz invariance 
of the Hamiltonian. It means that the considered nonperturbative approach for describing the 
internal structure of the "physical" electron should lead to the energy dispersion law (|43p for 
the entire range of the momentum P. 

The rigorous method of taking into account the translational symmetry in the strong coupling 
theory for the "polaron" problem was elaborated in the works of Bogoliubov [14] and Gross [19] . 
Let us remind that this method was based on the introduction of the collective variable R 
conjugated to the total momentum operator P. The canonical character of the transformation 
caused by three new variables Ri was provided by the same number of additional conditions 
imposed on the other variables of the system. In the "polaron" problem the quantum field 
interacting with the particle contributes to the total momentum of the system. It allows one 
to impose these conditions on the canonical field variables |14[ [T9] and the concrete form of 
the variable transformation is based mainly on the permutation relations for the boson field 
operators. 

The considered problem has some specific features in comparison with the "polaron" problem. 
Firstly, the formation of the one-particle wave packet is the multi-particle effect because this 
packet includes all initial states of electron-positron field as the fermion field. Secondly, its self- 
localization is provided by the polarization potential of the scalar field that does not contribute 
to the total momentum of the system. Therefore, we use a different approach in order to 
select the collective coordinate R. Let us return to the configuration representation in the 
Hamiltonian ([2]), where QED is considered to be the totality of N [N — ► oo) point electrons 
interacting with the quantum electromagnetic field in the Coulomb gauge [17] : 

N ~ t, i r - 

H = ^2{a a \p a + e A(f a )] + (3 a m + e <p{r a )} - - dr(V(p{r)) 2 + ^ uj(k)h- kX , 

kX 



a=l 

u(k) = k, % A = ctcn A , A = 1,2, 



"kX °fcA fcA' 

l(r) = V -7==e< A) [crs + c+ e"*] . (44) 

kX 

Here Q is the normalized volume; c ^ x ( c k\) are the operators of the creation (annihilation) of 

quanta of a transversal electromagnetic field, the quantum having the wave vector k, polariza- 
tion e( A ) and energy uj(k) = k. The sign of the interaction operators differs from the standard 
one because the parameter eo is introduced as a positive quantity. 
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As it was stated above the zero approximation of nonperturbative QED is defined only by 
a strong interaction of electrons with the scalar field, where the interaction with the transversal 
field is to be taken into account later in the framework of the standard perturbation theory. So, 
the conservation of the total momentum in the processes with transversal electromagnetic field 
will be provided automatically [2]. Therefore, while describing the quasi-particle excitation we 
should consider the conservation of the total momentum only for the system of electrons. In 
the considered representation it is defined by the sum of the momentum operators of individual 
particles 



P 



N 

E 

0=1 



p a , Pl*l(p)> = Pl*l(p)>. 



It means that in the configuration space the variable R conjugated to the total momentum 
is simply a coordinate of the center of mass of the electron-positron system and the desired 
transformation to new variables is as follows: 



f a = R + p a , 

P = -iV R , 
The Hamiltonian 



R 



1 

N 



N 



N 



Pa 



-*v« 



o=l 



0=1 



N 



VR+p' a , 



P'a 



N 



N 



6=1 a=l 

with new variables has of the following form 



H 



N 

E 



N 



V R + p' a + e A(p a + R) 



+ p a m + e o 0(p a 



- / dr(V(p{r))' 



+ uj(k)h 
k\ 



k\- 



(45) 



It should be noted that the matrix elements of an arbitrary operator in a new configuration 
representation are to be calculated by integration over the coordinates both of the center of 
mass and the relative variables. Let us introduce a special notation for this norm: 

(46) 



«$i|M|$ 2 » = / dR\[dp a §X(R,{p a })M§ 2 (R,{p a }). 



The interaction between the "physical" electron and the transversal electromagnetic field will 
be taken into account by means of the perturbation theory (see below Section [5]) and only the 
scalar field is considered in the zeroth approximation. Let us denote by Hq that part of the 
operator (|45p which does not depend on the transversal electromagnetic field and describes the 
internal structure of the "physical" particles in the zero approximation. In fact, the operator Hq 
does not depend on the coordinate R because of its commutativity with the operator of the total 
momentum of the system of electrons. This also follows from the well known result [T7] that in 
the Coulomb gauge the scalar potential could be excluded from the Hamiltonian. As a result 
the operator Hq depends only on the vector differences (r a — i%) = (p a — Pb) and does not change 
with the simultaneous translation of all the coordinates. As a consequence, the eigenfunctions 
of the Hamiltonian Hq depend on the coordinate R in the same way as for a free particle: 

1 



$(R,{Pa}) 



e iPR \^(P,{p a })}, 



Hq 



Po(P) 



(2vr) 3 / 2 
N 

= E 

a=l 



+ p a m + eoif(pa) 



^ / df(V0(f)) 2 
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Further calculations consist in returning to the field representation by the variables p a in the 
limit (N — > oo) and in using the approximate trial wave packet |3>i(P)) similar to @ but with 
the coefficient functions depending on P. Thus, in the framework of non-perturbation QED the 
orthogonal and normalized set of states for the one-particle excitation of the electron-positron 
field is defined as follows: 

l*SV)> = jAsJjj™/ ^(^J + ^^OiO;^). r = R + p, (47) 

with the norm (|46p and the coefficient functions Uff s (P), Vg s (P). 
Using the coordinate representation for these functions 

Mr,P) = J T^^E^K-^, = J (48) 

one can find the following functional for calculating the value E{P) corresponding to the energy 
of one-particle excitation with an arbitrary momentum 

E(P) = j dr^ + (f,P) (aP - iaV + /3m ) + e ^(f, P) #(r,P) 

- y +c (r,P)[(-aP-iaV + (3m ) + e ^(f, P)] }^ c (r, P), 

r dr' - -> 

(p(f,P) = e / — — — \^ + (f f ,P)^(f",P) - y +c (f",P)^ c (r",P)], 
J \r — r'\ 

J dr[^ + (r,P)^(r,P) + y +c (r',P)V c {r',P)] = 1. (49) 

It should be noted that the mean value of the operator Hq(P) in the representation of the 
second quantization was calculated by taking into account the expression 

hm Iy = i, 

p 

that follows from a well known formula for the density of states when EPF is quantized within 
the normalized volume f2 [2]. 

It was mentioned above, that there is quite a close analogy between the nonperturbative 
description of QED and the theory of strong coupling in the "polaron" problem. Therefore, it 
is of interest to compare the obtained functional (|49|) with the results of various methods of 
including the translational motion in the "polaron" problem. The simplest one was used by 
Landau and Pekar [35] who introduced the Lagrange multipliers in the form 

J(P) = J(P = 0) + (FV), (50) 

with the functional J(P = 0) referring to a resting "polaron" and the Lagrange multiplier Vj, 
denoting the components of the quasi-particle average velocity. 

We can see that the obtained functional (|49p has the same form as (j50j) if the relativistic 
velocity of the excitation is determined by the formula 



V = j dr{ty + (r,P)(a)^(r,P) + ^ +c (r,P)(a)^ c (r,P)}, 



corresponding to the well known interpretation of Dirac matrices [2]. 
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Varying this functional by taking into account the normalized conditions leads to the following 
equations for the wave functions 

{(aP - iaV + (3m ) + e ip(r, P) - E(P)}^(r, P) = 0, 

{{aP + iaV - f3m ) - e <p(r, P) - E(P)}^ c (f, P) = 0. (51) 

These equations show that unlike the "polaron" problem [TJ], the translational motion of 
the quasi-particle determining the momentum P is related in our case to its internal movement 
described by the coordinate f by means of spinor variables only. The physical reason for this 
separation of variables is explained by the fact that in QED the self-localized state is formed by 
the scalar field, and its interaction with the particle does not involve the momentum exchange. In 
order to find the analytical energy spectrum E{P) the system of non-linear equations (|5ip should 
be diagonalized with respect to to the spinor variables. The possibility of such diagonalization 
seems to be a non-trivial requirement for the nonperturbative QED under consideration. 

The solution of the equations (|51|) can be found on the basis of the states for which the 
dependence on the vector q in the wave packet amplitudes U^ s , V^ s remains the same as it was in 
the motionless "physical" electron. However, the relation between the spinor components of these 
functions can be changed. But as the self-consistent scalar potential involves the summation over 
all spinor components it does not depend on the momentum for the class of states in question: 

V>(r f ,P) = <p(r)\p =Q . (52) 

In the coordinate representation the transformation of the spinor components of the wave 
functions satisfying the equations (|5ip takes place because of the dependence on the momentum. 
The solution can be constructed by sorting out various linear combinations of the wave functions 
\j/(r), i|/ c (r), ^(r), \j/ c (r) found in Section [2] for a resting electron. These functions correspond 
to the degenerated states in the case of P = but are mixed for a moving electron. It is found 
that there is only one normalized linear combination satisfying all the necessary conditions of 
self-consistency: 

*(r,P) = L(P)¥(r) +K(P)^ c (f), ^ c (r,P) = Li(P)* c (r) + #i(P)#(f), 

|L| 2 + \K\ 2 = |Li| 2 + \K X \ 2 = 1. (53) 

The condition (|52j) . according to which the potential does not depend on the momentum for 
the excitation at the energy E(P), is fulfilled if the coefficients are related as 

L x = -K, K\ = L. 

These relations are also consistent with equations (|51j) for wave functions. 

This means that the "physical" electron moves in such a way that its states are transformed in 
the phase space of the orthogonal wave functions (|16p - (l24h but the amplitudes of its "internal" 
charge distributions are not changed. The analogous approach is considered for Lorentz invariant 
transformation of the bispinors corresponding to the free particles in QED |37j . 

Substituting the superpositions ([53]) into equations (|5T|) we use the following relations: 

(aP>= ( ^Xf^iAM' \ (aP)y c - ( 9i(r)(aP)Sl 1/2il)M \ 



-iaV + 



g(r)(aP)n 1/2fi ,M J ' \ -ifi(r)(aP)n 1/2i0! M 

V */(n n i/2,i,M / 

V 5n r )^i/2,M/ / 

and similar formulas for the functions ^(r), x I /C (r), <7j are the Pauli matrices. 
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For equations (|5ip to be fulfilled for any vector r it is necessary to set the coefficients of 
spherical spinors equal to the same indexes /. The corresponding radial functions are proved to 
be the same under these conditions, and the following system of equations for the spinors x^i 
is obtained 

iL{aP)xt + K(E + E )xf = 0, iK{BP)xt + L(E - E ) X t = 0, 

L(aP)xt + iK{E + E ) X ^ = 0, K(aP) X f + %L{E - E ) X f = 0, 

Li(?P)xI ~ iKi(E - E )xi = 0, K 1 (aP) X o - iL x {E + E ) Xo = 0, 

iL 1 (aP)xi -K^E- E )xi = 0, iK 1 (aP) Xo - L X {E + £ )Xo = 0, (54) 

where the notation E(0) = Eq was used. 

Spin variables in equations (|54p are also separated. In order to show this one can use, for 
example, the relation between the coefficients resulting from the 4th equation in (|54p in the 
first one of these equations: 

+ . K(aP)xt 
Xi = % - 



L(E - E 



As a result there exists a non-trivial solution of these equations for two branches of the energy 
spectrum 



E e>p = ±^jE$ + P 2 , (55) 

referring to the electron and positron zones, respectively [29]. The same expressions can be 
obtained for all conjugated pairs of the equations in (|54|) . The coefficients in the wave func- 
tions (I53j) can be found taking into consideration the normalization condition: 



L e =K\ = — P K e = -L\ Ee E ° 



VP 2 + We ~ Eo) 2 ^P* + (E e -E ) 2 

If = Kf= P =, KP = —L\ = Ee + E ° = , (56) 

y/P 2 + (E e + E ) 2 y/W+iET+W 

One can see that this set of coefficients coincides with the set of spinor components for 
solving the Dirac equation for a free electron with the observed mass m = Eq [37] • Thus, the 
results of this section show that the "internal" structure of the "physical" electron (positron) 
considered in this paper is consistent with the energy dispersion (|55p for a real free particle due 
to the relativistic invariance of the Dirac equation. 

For the interpretation of the wave packet (|47p as the state vector for a "physical" particle it 
is essential that it is the eigenvector of the total momentum of the electron-positron field in the 
framework of the considered zeroth approximation. It means that the "physical" electrons with 
different momenta form the set of orthogonal and normalized functions if the condition (|46p is 
taken into account: 

((^ 1 °\P 1 )\<S>? ) (P))) = 5(P-P 1 ). 

5 Perturbation theory for QED 

with the "physical" electron-positron field. Discussion 



It is well known that the renormalizability is one of the most important features of QED and it 
is confirmed by the coincidence of its results with the experimental data. Actually it means that 
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the calculated characteristics of the real electromagnetic processes do not depend on the initial 
parameters eo, tuq but only on the observed values of e and m. So, we should clearly show that 
the effects of the nonperturbative QED refer to the internal structure of the "physical" electron 
(positron) only. However, its interaction with real electromagnetic field should be defined by 
the standard perturbation theory with the fine structure constant a <^ 1. 

Besides, it should be demonstrated that the integrals which define the corrections to the zero 
approximation in the perturbation theory with charge e ~ 6q are converged without introducing 
any additional regularizating parameters including the cut-off momentum. 

To solve these problems, it is necessary to consider the form of the perturbation theory which 
uses the basis of states corresponding to the "physical" electron (positron) with various momen- 
tum. Taking into account the results of the last section let us represent the QED Hamiltonian 
([2]) in the following form: 

H = H + Hi, 

H = I dr\4*(^[-id(V R + V r ) + Pm ]^:+e ip(^:p(^ 



kX 

H I = e J dr:^*{f)[aA{r + R)]^{r):. 



The fact that the Hamiltonian depends on the coordinate R canonically conjugated to the 
total momentum of the electron-positron field determines the spontaneous breaking of symmetry 
of the system. This situation was discussed by Bogoliubov [27] for a lot of concrete physical 
systems. However, the global symmetry of the system is reconstructed when the observed values 
are averaged by the coordinate R in accordance with the norm (06]). 

In the previous sections the spectrum of the one-particle excitations of the electron-positron 
field was calculated approximately in the leading order in a power series of Qq . Let us now 
consider the matrix elements of the perturbation operator Hj that are defined by the transitions 
between the states corresponding to the "physical" electrons (positrons) with the momenta P 
and Pi and with emission (absorption) a quanta of the transversal electromagnetic field. It is 
important that the vacuum state of the "physical" electron-field is also renormalized with respect 
to to the "primary" vacuum of the "bare" particles. Actually this new vacuum is the wave 
function of the collective quasi-particle states localized in the vicinity of various coordinates R. 

Let us write the obvious form of the initial \i) and final |/) states of this process using the 
definition (I47p for the one-particle wave packet and the state vector of the free electromagnetic 
field corresponding to N quanta with the wave vector k and polarization A: 



^ I dq{U?A°)4s + VfrWb+yeW-^ ^ \0;0;<p(r- R);<p(f- Ri)), 



(27r) 3/ 2 y q ' s " 1 * s q ' s " ' ^ [-® 

kX 



\f) = 1^(0), *i(A)Jvg) = j^e iPlMl J dq{U<tAPiK + VfrlPOb+yeW-** 



- (cT ) 

j dq{U^ s {Q)a+ s + Vk(0)6+K*P-*> k >* |Q; 0; <p{f - R); <p{f - &)). 

k,X' 
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Let us calculate the matrix element of the operator Hj for the transition of the "physical" 
electron between the states with 4-momenta P = (P, E) and P\ = (Pi, -Ei) with the emission of 
one quantum of the transversal electromagnetic field with k = (k,uj). It is supposed that both 
fermion and photon 4-momenta are out of the mass surface. In this case the matrix element 
Mfi = T(P, k)5(P — P\ — k) defines the vertex function T(P, k) in the considered representation. 
The norm (1461) is used for calculation: 



is: 



Mfi = ((f\Hi\i)) = -|=EE/^/^/ d ^ iPn - Pvlil) ^ 

s,s' 

-^v(A)^,.'(o)^(o)^(P)] 

If the coordinate representation (|4'8|) is used for the coefficient functions U$ iS , V^ s , the result 

M H = -|== JdrJdHj dRi{[$>*(r- R, 0)(ae A )*(f - R, P)^*(r - R u A)*(r - Ri,0) 
+ ^*(r- R, 0)*(r - R, P)**(r - R x , Pi)(ae A )*(f - R x , 0)] 

- [f c '*(r- J R,0)(ae A )^ c (r-^,P)^ c '*(r- ^1,^)^(^-^1,0) (57) 
+ ^ c '*(r - P, 0)^ c (r-P, P)* c '*(r-i?i, Pl)(ae A )^ c (r-Pi, 0)] } e «P-R-PvRi) e -ikr 

Let us pass to the integration over the coordinate of the center of mass R c = {R + R\) /2 and 
relative coordinate Rq = (R — R{) in equation ([57]) and introduce the following vectors: 

ri = r - R c - -Ro, r 2 = r-R c + ^R . 
Then the vertex function results in: 

Mf i = T{P,k)S(P-P 1 -k), P x = P-k, 

r(P,k) = -JU J df 2 e ip ^ J df 1 e- < ^{**(r 2j P* 1 )(^)*(r 2 ,0)**(7' 1> 0)*(f > 1) ^) 

+ ^* (r 2 , Pi ) (f 2 , 0) * (n , 0) (ae A ) ^ (n , P) 
-^*(r 2 ,Pi)(ae A )^ c (r 2 ,0)* c '*(fi,0)^ c (fi,P) 

- * c '*(f 2 ,A)* c (r 2) 0)^*(n,0)(ae A )* c (ri,P)}. (58) 

The linear combinations (j53[) and the wave functions ^ ^f c from equations ()16p . (I20p . 

(|24p . (|25p for the resting "physical" electrons should be used in equation (|58p . Then one can 
pass to the two-component spinors in the integrand and find the following formula: 

T = i-^=j dn j ^{p(rl)p(r 2 )( X +*(d ; e A )xo + )[^*(A,Pi)MP^) 

-L*(Pi,Pi)K(P,P)]-pi(ri) /0 i(r 2 )( X o*( ( r e -* A )xo") 

[Kt{P ll E 1 )L 1 {P,E) - L* 1 (P 1 ,E 1 )K 1 (P,E)]}e^e-^, 



x 



where the functions L(P,E), K(P,E), L X (P,E), K\{P,E) are defined by equations ([56]) but 
with the energy variable out of the mass surface E ^ Eq. The charge densities are defined as 
follows: 

p{r) = g\r) + / 2 (r), Pl {r) = gf(r) + ff(r). 
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Accordingly to Section [2] the functions /, g and f±, g\, that describe the "physical" electron 
internal structure, can be expressed over the universal function and are differed from each other 
only by the normalization conditions 



1 3/2«o(™7i ) , > C 3/2Uo(rm ) 

TTc m ° — T" ' 9l{r) = TTc m ° 



\ 1 3/2V (rm ) C 3/2V (rm ) 

f{r) = TTc m ° ~ ^~ ' 9l{r) = TTc m ° — T- 

Besides, the spin and angular variables can be excluded by means of equations (I53j) for the 
spinor functions of the "physical" electron. As for example, 

Taking into account the normalization conditions for the spinors Xo~'~ one can find: 

(Pe A ) (Pie A ) 



E + £7(0) Ei + £7(0) 



F*(|P-fc|)F(|P|)], 



F(\P\) = F(a, \P\/m) = / ds[t^(x) + i^(x)] 



2 , sin \P\x/mo 



o |P|x/mo 



oo 



2 , sin a|P|x/2|ao|?n 



dx[^(x)+^(x)] ' ' ^ 1 Ul ■ (59) 
o a\P\x / 2\ao\m 

If one use the relations eo(l — C)/(l + C) = e, (1 — C)/(l + C) ~ e^" 2 , £7(0) = m and put 
£(|P|) = 1 the expression (|59p will coincide completely with the analogous matrix element in 
standard form of QED with the renormalized charge e and mass m of the electron [17] . It means 
that the series for the observed characteristics of the "physical" electron in a reverse power of the 
"primary" charge e$ 1 leads to the terms with the same structure as the standard perturbation 
theory in a power of e. 

However, equation (|59p shows that behavior of the renormalized vertex function r(P, k) at 
large arguments is defined by the Fourier transformation of the radial functions that describe 
the charge distribution in the "physical" electron. Analysis of the asymptotic behavior of the 
solutions of equations (128j) shows that the form- factor P(|P|) at large |P| S> is described as 
follows: 

F(\P\) =r-^ , L = 2?n = m-^, (60) 

U U (|P| 2 + £2)2' a K ' 

and tends to zero rather quickly. 

So, if the representation with the basis of states corresponding to the "physical" electrons was 
used in the zeroth approximation then the cut-off momentum Lq ~ 1 GeV appears naturally 
and it has definite physical meaning. It means that in this representation the perturbation 
theory in a power of the observed value of the fine structure constant a< 1 does not contain 
the ultraviolet divergencies in the integrals. 

It may seem at first sight that the finite value of the cut-off momentum is in the contradiction 
with the experiment. It is well known that the comparison of the theoretical calculations in the 
framework of the standard renormalized QED with the experimental data even in the low energy 
range (for example, for anomalous magnetic moment) shows that no additional momentum 
L < 10 7 GeV can be introduced in QED (for example, |33j). However, these results do not 
maintain that the electron is the point-like object with the form- factor equal to unity in the 
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entire range of the momentum because the form-factor of the electron appears due to the high 
order corrections of the perturbation theory with any calculation scheme [2J. Actually it is 
followed from these results that the electron cannot be considered as the composite system 
including any parameter L, as for example, the hypothetical particle with a mass M ~ L [33] . 

But the most essential feature of the present approach is the result that the cut-off momen- 
tum Lq is not the additional parameter of the theory because it can be expressed only in terms 
of the observed values of QED, namely a and m(\). The form factor of the "physical" electron 
in this case can be considered as the result of the partial summation of high-order corrections 
of usual perturbation series. Let us now show that the appearance of Lq at the intermediate 
calculation stages may not contradict to the fact that it is not effect at the final results for the 
observed physical values in the renormalised QED. 

Really, let us consider a general structure of the non-renormalized form factor of the electron 
as the function of the invariant parameters in the framework of QED with an arbitrary cut-off 
momentum L and the non-renormalized values of ao, niQ. For definiteness, let us write only its 
scalar part: 



Jv^ 7TZ-Q 

F (k) = F ( a ; — 



(61) 



Computation of this form factor by means of the perturbation theory with the standard 
renormalization scheme is defined by the calculation of the vertex function [33j and corresponds 
to the expansion of the function (|6ip as a power series in both coupling constant ao an d the 
parameter ijiq/L: 



F (k) » l + a 



? (o) 

01 



F, 
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02 



0; 



k 2 
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+ L *01 



k 2 
0;— ;0 
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02 
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k 2 
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2 ' 
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k 2 
0;— ;0 



/?)? 
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+ 



(62) 



Let us remind that the QED renormalizability means that all coefficients of these series are 
remained finite in the limit L — ► oo if the parameters ao, mo are changed by their renormal- 
ized values a, m. Then the theoretically calculated value of the form factor is defined by the 
expansion in a power series of both a and an independent parameter £ ~ m/L: 



F th (k) « 1 + aF^ 



(0) 
01 



m 2 



+ cfF. 



2£,(0) 



02 



k 2 



m- 



+ 



m 
1 



aF, 



(i) 

01 



m 2 



+ « 2 4 1} ( 



(k 2 



rn- 



+ 



+ 



(63) 



The formal series that appear in the intermediate stages and define the relations between the 
"primary" parameters of the theory ao, tuq and the observed value a, m can be represented as 
follows: 



Q 



m 



ao, ^) , 



m 



1 + ao^x 



+ a$B 2 



m \ 

a Zi [ a , — J 



(64) 



Actually all coefficients of these series tend to infinity in the limit L — > oo because of the 
divergent integrals. Therefore the concrete values ao, tuq are indefinite if usual renormalization 
scheme is used. 

Analysis of [53] shows that in the result of the renormalization the additional terms contained 
the parameter £ in the expansion ([63]) are absent at least with an accuracy ~ 10~ 10 corresponding 
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to L < 10 7 GeV. Analogous estimation is correct when calculating characteristics of other 
electromagnetic processes in QED. It means that after passing to the limit L — > oo the formulas 
of type of equation (|63|) give the practical algorithm for calculation of any observed values in 
QED with very high accuracy. For that purpose the internal inconsistency of the theory that is 
conditioned by indetermination of the series (|64p becomes unessential. 

We believe that the variational approach presented in the paper allows one to preserve the 
computing possibilities of QED but to avoid its logical inconsistency. As it was shown above, in 
this approach the "physical" electron characteristics are defined by the quasi-particle collective 
excitation of the "primary" electron-positron field itself and are not conditioned by additional 
parameters or bond with any hypothetical particle. Therefore the perturbation series that 
define the corrections to the form factor conditioned by the interaction with the transverse 
electromagnetic field are differed from equation ([62]) qualitatively According to equation ([59]) 
for the vertex function in this case the expansion parameter is the value Oq 1 ~ a and the form 
factor depends on the invariant parameters as follows: 



k 2 



m 



F(k) = F[a ;—;—\ 
m z mo / 



a 



a 



mo ~ m/a, 



and the effective cut-off momentum for the divergent integrals is the value Lq ~ mo 
Then the following series arises instead of the expansion ([62]): 
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(65) 



In order to calculate the observed value for the form factor equation (|65p should be completed 
by the series that express the "primary" parameter atQ, mo over the observed values a, m. These 
series have the following general form: 



a « a 1 [Ao + a X X\ + 
m w m aQ 1 [[iQ + a 



i A 2 + ---" 



a Z(ao), 



o Vi + "o 2 / i 2 H ] = moZ^ao). 



(66) 



According to equation (|60|) all coefficients in these series are finite and and do not depend on 
additional parameters. 

Substitution of equation ([66]) in equation ([65]) shows that the electron form factor calculated 
by means of the considered approach contains only the finite values and does not depend on the 
additional parameters as it should be in accordance with the analysis in [33J: 
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The considered qualitative analysis of the "physical" electron form factor can be applied not 
only for the interaction with the vacuum transversal electromagnetic field but also for any kind 
of the external fields. It means, for example, that after renormalization with the self-localized 
electron states the QED corrections for the "physical" electron bound in the Coulomb field of 
the nucleus in the hydrogen-like atom should be described by the same formulas that can be 
found with usual renormalization scheme |36|. 
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Actually it means that the considered form of QED that uses the basis of the self-localized 
states and the perturbation theory on the parameter 1 can lead to the same results for the 
observed values as the renormalized series in the physical fine structure constant a. They are 
differed only when calculating the connections between the primary and observed values of 
«o an d ol as well as mo and to. However, the strict proof of the quantitative coincidence of 
our results with known ones are connected with rather cumbersome calculations of high-order 
corrections in a power of the parameter 1 on the basis of the vertex function (|59p and the 
propagation function for the "physical" electron and should be considered separately. 

In any case it is evident that for the practical purposes the previously existing calculation 
scheme in QED is more efficient than one considered in the present paper. Therefore the proof 
of the equivalence of the both forms of QED is important mostly from the point of the internal 
self-consistency of QED. Nevertheless, the suggested approach may have the applied interest as 
well because it allows one to consider other leptons as excitations of the "primary" electron- 
positron field (see Section [3] for ^-meson). It is also may be essential for the constructive use of 
the renormalization group in QED and for the nonperturbative analysis of other quantum field 
models including non-renormalized ones. 

6 Conclusions 

Thus, in the present paper the self-localized quasi-particle excitation of the electron-positron 
field is found for the first time in the framework of a standard form of the quantum electro- 
dynamics. Existence of a new state in the spectrum of the well investigated Hamiltonian is of 
great theoretical interest by itself. 

Besides, the physical interpretation of the obtained solution as the "physical" electron has 
been considered. It allowed one to find in the finite form the relation between the charge eo 
and the mass iuq of the "bare" electron considered as the "primary" parameters of the theory 
and observed values of the charge eo and the mass too of the "physical" electron. The obtained 
relations between these values (eo ~ e , mo ~ me -2 ) show that they cannot be calculated by 
means of the standard perturbation theory. 

The important consequence of the considered interpretation of the self-localized state is the 
possibility to represent the Dirac idea and to consider /i-meson as an additional quasi-particle 
excitation of the electron-positron field. It is shown that the calculated mass of //-meson is very 
closed to its experimental value. 

It is also shown that the self-localized quasi-particle state does not contradict to the Lorentz- 
invariance of the system and its binding energy can be considered as the observed mass of 
this excitation. It permitted us to describe the interaction of the "physical" electron with 
the transverse electromagnetic field and to show that it can considered on the basis of the 
perturbation theory as the series in a power of the observed charge of electron e ~ e^ 1 <C 1. It 
proved that the coefficients of these series are free from the ultraviolet divergence. 

The internal structure of the "physical" electron considered as the collective excitation does 
not depend on any additional parameters with the exception of a and to. It is described 
by the finite functions without introduction of the cut-off momentum L. It is also discussed 
quantitatively that an existence of such structure does not contradict to the experimental data 
about the electron form factor. 
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